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THE  PLASMA  ASSISTED  MODIFIED  BETATRON 

I.  Introduction 

In  order  to  obtain  a  high  current,  high  energy  electron  beam,  several 

1  -3 

laboratories  are  giving  serious  consideration  to  the  modified  betatron. 

The  idea  is  that  the  net  space  charge  forces,  which  is  outward  in  the  poloidal 
plane,  can  be  confined  with  a  strong  toroidal  field.  Theory  and  particle 
simulations  show  that  for  at  least  one  poloidal  drift  orbit,  much  higher 
currents  can  be  confined  in  a  modified  betatron  than  in  a  conventional 
betatron.  Self-consistent  fluid  formulations* have  shown  that  cold  fluid 
equilibriums  exist  and  can  be  accelerated.  Recently,  the  theory  was  extended 
to  include  transverse  pressure.^ 

However,  while  the  modified  betatron  offers  many  advantages,  it  is 
possible,  but  not  certain,  that  future  experiments  may  face  several 
difficulties.  First  of  all,  one  must  Inject  the  beam  across  toroidal  field 
lines  in  order  to  have  a  beam  centered  in  the  liner.  The  current  scheme7 
proposes  to  shoot  the  beam  into  the  toroidal  vacuum  chamber  near  the  liner. 

The  drift  due  to  the  focusing  fields  and  image  fields  cause  the  beam  to  drift 
in  the  poloidal  plane  around  the  liner.  In  one  toroidal  transit  (about  20 
nsec)  it  should  drift  enough  to  miss  the  injector.  In  one  poloidal  drift  time 
(several  hundred  nanoseconds),  external  fields  can  be  changed  to  bring  the 
beam  slightly  in  from  the  liner  so  that  it  misses  the  injector  again  and 
henceforth.  On  a  longer  time  scale,  wall  resistivity  causes  the  beam  to  drift 
inward. 

However  there  is  a  significant  range  of  beam  currents  for  which  wall 

resistivity  causes  the  beam  to  drift  outward  if  it  is  near  the  liner,  but 

Inward  if  it  is  near  the  center.  Since  it  is  unlikely  that  the  beam  can 

reverse  drift  directions  on  the  way  in,  the  injection  scheme  of  Ref.  7  appears 

to  be  viable  only  for  fairly  low  beam  currents. 

Manuscript  approved  October  18,  1984. 


1 


A  second  possible  difficulty  concerning  the  modified  betatron  is  the 

diamagnetic  to  paramagnetic  transition. Depending  on  whether  the  net  self 

force  is  outward  or  inward  in  the  poloidal  plane,  the  net  electron  drift 

velocity  is  in  the  diamagnetic  or  paramagnetic  direction.  As  the  high  current 

beam  accelerates,  it  makes  a  transition  from  diamagnetic  to  paramagnetic 

current  flow.  It  has  been  shown  that  subject  only  to  the  constraint  that  the 

-3 

acceleration  time  ~  10  sec  is  very  long  compared  to  the  drift  time,  > 
10"  sec,  this  transition  must  suddenly  change  the  topology  of  the  beam  orbits 
in  the  poloidal  plane. Whether  the  beam  can  survive  such  a  sudden,  violent 
perturbation  is  an  open  question. 

Finally,  although  the  focusing  fields  in  the  modified  betatron  stabilize 
the  l  *  1  resistive  wall  instability,  l  m  2  modes  are  still  unstable  and  pose 

O 

a  real  threat  to  beam  confinement  in  the  modified  betatron.0  These 
difficulties  are  discussed  in  more  detail  in  Sec.  II. 

The  root  cause  of  all  of  these  problems  is  that  the  beam  self  fields  in 
the  modified  betatron  are  outward  in  the  poloidal  plane.  By  making  the  self 
fields  inward,  all  three  of  these  difficulties  can  be  alleviated.  One  way  of 
doing  this  is  with  a  very  low  density  background  plasma.  This  is  discussed  in 
Sec.  III. 

Section  IV  discusses  other  physics  issues  concerning  the  plasma  assisted 
modified  betatron,  among  them  plasma  production,  plasma  response  to  the 
injected  beam  and  ion  resonance  and  streaming  instabilities. 


Injection,  Diamagnetic  Transition  and  l  ■  2  Resistive  Wall  Instabillt 


in  a  Vacuum  Modified  Betatron 

This  section  discusses  what  appears  to  be  three  fundamental  Issues 
regarding  the  high  current  modified  betatron  operating  with  a  vacuum 
background:  beam  injection,  the  diamagnetic  to  paramagnetic  transition,  and 
the  t  m  2  resistive  wall  instability. 

A.  Injection 

One  of  the  important  issues  for  the  modified  betatron  Is  Injecting 
the  beam.  The  present  thinking  for  the  NRL  modified  betatron  experiment 
is  described  in  Kapetanakos,  et  al.^  The  beam  is  Injected  near  the  liner 
and  drifts  around  the  edge  of  the  liner  through  a  combination  of  drift 
generated  by  the  focusing  fields  (field  index)  and  image  forces.  The 
former  is  directed  Inward  in  the  poloidal  plane,  the  latter,  outward.  At 
high  beam  current  the  latter  dominates  and  the  beam  rotates  in  say  a 
counterclockwise  direction  as  in  Fig.  3  of  Ref.  7  for  the  case  of  a  ten 
kilo  Amp  beam.  If  the  combination  of  forces  is  large  enough,  the  drift 
velocity  will  be  great  enough  so  that  after  one  toroidal  revolution,  the 
beam  will  be  displaced  in  the  poloidal  plane  by  a  large  enough  distance 
that  it  misses  the  injector.  Then,  since  it  has  many  more  toroidal 
transits  before  it  would  hit  the  injector  again,  macroscopic  fields  could 
change  sufficiently  to  bring  the  beam  into  the  center. 

One  potential  problem  with  this  scheme,  which  is  not  addressed 
in  Ref.  7  is  that  for  high  current  beams,  the  net  poloidal  force  on  the 
beam  is  outward  near  the  liner,  but  inward  when  the  beam  is  at  the 
center.  Thus,  as  the  beam  continues  to  spiral  in  the  poloidal  plane,  at 
some  point  it  must  reverse  direction.  To  see  this  more  quantitatively, 
if  the  field  index  of  the  beam  is  4,  the  focusing  field  produces  an 
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inward  poloidal  force  on  a  charge  q  of  q  *■  where  RQ  is  the  major 


radius  of  the  equilibrium  orbit,  B_  is  the  vertical  field  and  p  is  the 


displacement  of  the  beam  from  the  equilibrium  orbit  in  the  poloidal  plane 


2  2  2 

o  *  (R  -  R  )  +  z  .  The  image  electric  force  for  a  cylindrical  system 


is  given  by 
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2l(Amps)  p(cm) 


10(a2  -  p2) 


where  a  is  the  minor  radius  of  the  liner. 


The  actual  image  force  is  canceled  in  part  by  magnetic  forces  and 


also  by  any  fractional  neutralization  f.  If  the  beam  is  near  the 


wall  (p  *  a)  the  net  poloidal  drift  velocity  is  given  by 
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C  KY  t}  106(cm)  b‘  (Gauss)  "  2Ba  R 
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where  6  is  the  distance  from  the  beam  center  to  the  liner.  Since  the  beam 


enters  the  toroidal  liner  right  near  the  outer  edge  of  the  liner,  i  is 


roughly  equal  to  the  beam  radius  p^. 


If  the  beam  is  near  the  center  (so  p  <  <  a),  the  poloidal  drift  is 


given  by 
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For  the  case  of  a  vacuum  modified  betatron,  the  current  can  be  classified  as 


being  in  one  of  three  ranges: 
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I.  High  Current 


2  B  a 

5  a(cm)Y  BQ  (Gauss)  1(A)  (4) 
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II.  Intermediate  Current 

2  B  a  -  B  a 

10 6(cm)BQ  (Gauss )y  (gy— g-)<  1(A)  <  5YVcm)B6  (Gauss)  fy  )  (5) 
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III.  Low  Current 

2  B  a 

1(A)  <  10y  <$(cm)  B0  (Gauss)  f2— z—  - -).  (6) 

^D0  o 

In  the  high  current  regime,  the  forces  on  the  beam  in  the  poloidal 
plane  are  outward  and  the  beam  always  rotates  in  the  counterclockwise 
direction.  In  the  Intermediate  regime,  the  forces  are  outward  when  the  beam 
is  near  the  wall,  but  inward  when  the  beam  is  near  the  center.  Thus  in  this 
current  regime,  the  beam  must  reverse  its  direction  of  rotation  before  it  gets 
to  the  center.  Also,  at  some  radius  between  the  center  and  the  wall,  the  beam 
will  have  zero  poloidal  drift  velocity.  It  seems  likely  that  some  time  after 
injection,  an  intermediate  current  beam  will  stagnate  around  this  point  and 
gradually  fill  the  chamber.  In  the  low  current  regime,  the  inward  focusing 
forces  always  dominate  and  the  beam  rotation  is  clockwise. 

It  Is  also  worth  noting  that  if  the  liner  is  resistive,  the  beam 
will  spiral  Inward  if  the  net  force  is  inward  and  visa  versa.  Thus  a 
resistive  wall  can  only  trap  a  class  III  low  current  beam.  It  is  possible 
that  an  intermediate  current  beam  can  be  trapped  if  it  can  be  brought 
sufficiently  near  the  center  that  the  net  forces  are  inward.  For  the 


2  x  103, 
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parameters  of  the  NRL  modified  betatron  experiment,  »  B.  * 

2 

Y  ■  4,  Rq  »  10  ,  a  *  15,  5*2,  *  140,  the  lower  and  upper  currents  of  the 

o  y 

intermediate  range  are  3.2  x  10 JA  and  1.2  x  104A.  Thus  the  maximum  current 
which  can  be  trapped  by  wall  resistivity  in  a  vacuum  modified  betatron  is 
about  3  kA.  Actually  however,  the  maximum  current  is  much  less  because  at  the 
3.2  kA  level  the  poloidal  drift  is  zero,  so  the  beam  will  strike  the  injector 
after  one  toroidal  transit.  Note  that  the  model  here  assumed  a  cylindrical 
system.  It  has  recently  been  shown*®  that  toroidal  effects  can  complicate 
this  picture  somewhat  by  allowing  more  complex  banana  shaped  orbits.  The 
authors  of  Ref.  10  suggest  an  injection  scheme  where  the  beam  orbit  is  shifted 
from  a  banana  to  circular  orbit  by  a  pulsed  change  in  toroidal  field. 

B.  Diamagnetic  to  Paramagnetic  Transition 

Let  us  imagine  that  the  beam  has  been  injected  and  is  centered  in 
the  modified  betatron.  The  question  then  is  the  individual  particle  orbits  in 
the  beam.  Each  particle  feels  an  inward  force  due  to  the  focusing  fields  and 
an  outward  force  due  to  the  self  fields.  If  the  latter  dominates,  the 
particle  has  an  F  x  B  drift  in  the  counterclockwise  direction,  analogous  to 
the  counterclockwise  whole  beam  drift  for  an  outward  image  force  discussed  in 
the  previous  section.  Then  the  J  (poloidal)  B  (toroidal)  force  is  inward.  In 
this  case,  the  beam  is  said  to  be  diamagnetic.  This  is  analogous  to  the 
terminology  in  plasma  physics,  where  the  poloidal  current  is  diamagnetic  if  it 
gives  an  Inward  force.  On  the  other  hand,  if  the  focusing  force  dominates, 
the  F  x  B  drift  is  clockwise  and  the  beam  is  said  to  be  paramagnetic. 

If  the  beam  has  uniform  density  and  radius  p^,  the  outward  force  is 
the  electrostatic  force  canceled  by  the  magnetic  force  and  fractional  charge 
neutralization.  A  test  charge  q  at  p  ■  p.  feels  an  outward  force 

I  O 
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The  inward  focusing  force  is  given  by 


p  Bz  °b 

f£  ’  - 11  -r 


so  Che  condidon  for  a  paramagnetic  beam  is  for  a  vacuum  modified  betatron  is 


-2  I  (AMPS)  ,  V°“SS>  % 

1  A  _  f _ \  '  -ws.  • 


10  rb(cm) 


Mote  that  a  high  current  beam  is  generally  diamagnetic.  However  as  it 
accelerates,  the  left  hand  side  becomes  smaller  as  y  increases,  and  che  right 
hand  side  becomes  larger  because  Bz  is  proportional  to  y.  Thus  for  a  high 
current  beam  which  starts  out  diamagnetic,  as  it  accelerates  it  ultimately 
makes  a  transition  and  becomes  paramagnetic. 

One  might  think  this  simply  means  that  the  poloidal  rotation  of  the 
particle  stops  and  changes  direction.  Actually  the  situation  is  considerably 
more  complex,  and  also  worse  from  the  point  of  view  of  operation  of  the 
modified  betatron.  In  a  recent  series  of  papers, it  has  been  shown  that 
subject  only  to  the  constraint  that  the  acceleration  time  is  very  long 
compared  to  the  drift  time,  an  approximation  well-satisfied  in  the  NRL 
modified  betatron  (but  not  satisfied  at  all  in  particle  simulations  of  the 
device),  the  diamagnetic  to  paramagnetic  transition  necessarily  results  in  a 
change  of  topology  of  the  beam.  The  outer  beam  particles  first  become 
paramagnetic  and  in  doing  so  scrape  off  the  edge  of  the  beam  and  form  a  large 
minor  radius  hollow  beamlet.  As  the  energy  continues  to  increase,  the  scrape 
off  point  moves  inside  the  beam  and  inner  beam  particles  continue  to  add  to 


the  outside  of  the  hollow  beamlet.  The  process  is  completed  when  the  beam  has 


turned  itself  completely  inside  out  and  has  gone  from  a  solid  to  a  hollow 
beam. 


Although  this  process  is  complicated,  it  is  verv  easy  to  see  that  in 
making  the  transition,  the  beam  must  turn  itself  inside  out.  To  show  this,  it 
is  only  necessary  to  invoke  the  conservation  of  toroidal  canonical  momentum 
Pq  .  If  the  poloidal  magnetic  field  is  given  by  Vt l>  x  i^  /R,  then 


PQ  -  YmRv0  +  (10) 

2 

To  evaluate  P0  ,  note  that  y  *  (E  -  q<J>)/mc  where  <J>  is  the  electrostatic 
potential..  For  a  cylindrical  beam  of  radius  r^, 


<fr(p) 


imqp 


P  <  P, 
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The  flux  has  three  components.  First  there  is  the  flux  of  the  vertical 
field  itself,  assumed  uniform 


B  R 
z 


(12) 


Secondly,  there  is  the  flux  associated  with  the  focusing  field.  If  the  field 
index  is  4,  this  is 


Note  that  ^  has  this  form  both  for  p  <  p^  and  p  >  p^.  Finally,  there  is  the 
flux  associated  with  the  self  field, 

•i>  -  (Va/c)2  R6(p).  (14) 

Thus  if  V  *  c,  near  the  axis  (p  *  0)  one  has  that 

2  2 
qB  p  nR  q  ir  , 

pe  -  -rf—  -  — ^ —  »  +  (15) 

Y  c 

Since  qBz  <  0  for  the  modified  betatron,  one  has  the  result  that  if  n  is  large 
enough  that  the  second  term  dominates  (that  is,  if  the  beam  is  diamagnetic), 

P.  (p  =  0)  is  a  relative  minimum.  However  far  from  the  beam  P„  is  dominated 
by  the  focusing  forces  which  have  the  opposite  sign.  Thus  PQ  as  a  function 
of  p  for  a  diamagnetic  beam  is  as  shown  in  Fig.  la.  On  the  other  hand,  if  the 
beam  is  paramagnetic,  the  first  term  on  the  right  hand  side  of  Eq.  (15) 
dominates  so  PQ  (p  *  0)  is  a  relative  maximum,  and  PQ (p)  is  shown  in  Fig.  lb. 

The  crucial  point  now  is  that  in  a  configuration  which  has  9 
symmetry,  P  is  an  exact  constant  of  motion.  Consider  then  the  orbits  at 
o  =*  0  and  p  =  pb  for  a  diamagnetic  beam.  The  former  is  inside  the  latter 
and  has  a  lower  value  of  P0  according  to  Fig.  la.  After  diamagnetic  to 
paramagnetic  transition,  these  values  of  P„  cannot  change.  However  a 
paramagnetic  beam  has  the  reference  orbit  at  a  relative  maximum  so  that  this 
must  correspond  to  the  orbit  initially  at  p  =  p^  in  the  diamagnetic  beam. 

Thus  in  making  transition  the  beam  must,  at  the  very  least,  turn  itself  inside 
out. 

Actually,  as  shown  in  Refs.  4-6,  not  only  does  the  beam  turn  itself 
inside  out,  it  transitions  from  a  solid  to  hollow  beam.  In  doing  so,  the  beam 
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could  strike  the  wall  and  thereby  disrupt.  Conditions  for  the  beam  to  remain 
confined  on  transition  are  given  in  Refs.  4-6.  However  even  if  the  beam  does 
remain  initially  confined,  it  is  not  certain  it  can  remain  confined  long  after 
suffering  such  a  violent  perturbation.  If  nothing  else,  the  hollow  profile 
produced  is  diocotron  unstable. 


Finally,  we  point  out  that  particle  simulations  of  the  beam  passing 
through  transition  have  been  done.  These  are  discussed  in  the  Appendix. 


C.  Resistive  Wall  Instability 

One  other  potential  difficulty  with  the  modified  betatron  is  the 
resistive  wall  instability.  If  a  beam  of  density  n  and  radius  o^  centered  in 
a  cylindrical  tube  of  radius  a,  the  frequency  of  a  perturbation  at  frequency 

O 

varying  like  exp  i2$  is 


oi  =* 


f(i  “  1)  +  f~)2* 

a 


(16) 


where  V  is  the  rotation  frequency  of  the  electrons 


V  -  ~  (y2  ~  f)  2imqp 

*  Be 


(17) 


and  oi^  is  the  frequency  of  rotation  generated  by  the  focusing  fields, 
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The  focusing  fields  produce  a  rotation  In  the  negative  9  direction.  Since 
q  <  0,  the  rotation  of  the  beam  itself  is  in  the  positive  9  direction  for  the 
case  of  a  vacuum  beam,  f  ■  0.  The  sign  of  the  frequency  is  such  that  as  long 
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as  a)  >  0,  wall  resistivity  gives  rise  to  growth  of  this  mode.  This  can  be 
understood  by  noting  that  since  the  beam  has  only  negative  charge,  the  net 
force  from  any  perturbation  must  be  outward.  Thus,  as  mentioned  in  Sec.  Ill 
A,  wall  resistivity  will  cause  the  beam  to  spiral  outward,  corresponding  to 
instability.  Since  the  natural  frequency  of  the  l  m  1  mode  is  very  low,  the 
sign  of  this  frequency  can  be  changed  by  the  focusing  fields,  thereby 
stabilizing  this  mode.  The  condition  is  that  the  beam  current,  as  defined  in 
Sec.  II  A  be  in  the  low  or  intermediate  regime.  However  the  t  ■  2  mode  has  a 
significantly  larger  frequency  so  that  in  the  vacuum  modified  betatron,  it 
cannot  be  stabilized  by  the  focusing  fields. 


*  «  *  .  '  •  n  «  •  «  •  ,  1  . 
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III.  A  Plasma  Assisted  Modified  Betatron 


In  the  previous  section  we  have  seen  that  the  modified  betatron  is  beset 

by  three  potentially  very  serious  problems:  the  difficulty  of  injection  at 

high  or  intermediate  range  currents,  the  problem  of  the  diamagnetic  to 

paramagnetic  transition  which  occurs  even  well  into  the  low  current  range,  and 

the  problem  of  l  >_  2  resistive  wall  instability.  However  all  of  these 

-2 

difficulties  resulted  from  the  fact  that  y  -  f  >  0,  or  that  the  beam  self 

forces  are  outward.  One  possible  cure  for  all  these  problems  then  is  to 

operate  the  high  current  modified  betatron  in  the  presence  of  a  low  density 

-2 

background  plasma,  so  that  y  -  f  changes  sign.  This  changes  the  sign  of  the 

self  forces  in  the  poloidal  plane  from  outward  to  inward  and  provides  cohesion 

for  the  beam  itself.  On  injection,  the  image  forces  are  now  toward  the  center 

so  that  wall  resistivity  can  trap  the  beam  in  the  center.  The  beam  will 

always  be  paramagnetic,  so  the  problems  of  making  transition  are  eliminated. 

Also,  the  frequency  of  the  diocotron  mode  will  be  negative  so  that  the 

resistive  wall  instability  will  be  stabilized.  The  plasma  densities  required 

are  low.  For  Instance  for  a  10  kA,  2  cm  radius  beam  with  y  ■  4,  background 

10  —3 

ion  density  of  order  10iU  cm  J  is  required.  For  a  1  kA  beam,  it  is  an  order 
of  magnitude  lower. 

Although  we  envision  the  betatron  as  operating  in  the  presence  of  a 

background  plasma,  the  scheme  proposed  here  has  little  in  common  with  a 

"plasma  betatron".  There,  the  beam  is  formed  from  runaway  electrons  and  the 

11—13 

beam  density  is  small  compared  to  the  background  plasma  density.  An 

alternate  scheme  involves  injection  of  an  astron  gun  produced  electron  ring 
into  a  high  density  collisional  plasma.^  In  the  scheme  proposed  here,  a  beam 
is  still  externally  injected,  and  the  preformed  plasma  density  is  low  compared 
to  the  beam  density. 
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IV.  Physics  Issues  Relating  to  the  Plasma  Assisted  Modified  Betatron 

In  this  section  we  consider  three  physics  issues  relating  to  the  proposed 
plasma  assisted  modified  betatron.  They  are  producing  the  low  density  fully 
ionized  plasma  in  the  toroidal  system,  the  response  of  this  plasma  to  the 
injected  beam,  and  the  ion  resonance  and  streaming  instability. 

A.  Producing  the  Background  Plasma 

10  -3 

Producing  a  fully  ionized  plasma  at  a  density  as  low  as  101U  cm  J 
appears  to  present  some  experimental  difficulties.  The  most  likely  approach 
would  be  to  produce  the  plasma  at  much  higher  density  and  let  it  expand  into 
the  entire  toroidal  chamber.  The  fact  that  there  is  a  vertical  field  which  is 
typically  five  to  ten  percent  of  the  toroidal  field  means  that  the  interior  of 
the  chamber  is  accessible  along  a  field  line  from  the  outside.  One  possible 
scheme  would  then  be  to  distribute  a  large  number  of  very  small  plasma  guns 
along  the  top  of  the  liner.  These  could  be  made  to  fire  simultaneously  so 
that  plasma  would  line  the  top  of  the  chamber.  As  it  drifted  in  along  the 
field  the  density  would  decrease  due  to  the  expansion.  The  plasma  would 
expand  along  the  field,  (and  also  outward  in  major  radius)  filling  the 
torus.  Since  the  expansion  velocity  is  about  10^  cm/sec,  it  would  take  many 
microseconds  for  the  torus  to  fill.  However  this  is  a  very  long  time  compared 
to  the  20  nsec  transit  time  of  the  beam  in  major  radius.  Thus  when  plasma 
conditions  are  optimum,  the  beam  would  be  fired  in.  The  best  way  to  determine 
the  optimum  plasma  configuration  is  almost  certainly  experimentally,  with  a 
small  scale  experimental  program. 


B 


The  Plasma  Response  to  the  Injected  Beam 

The  system  envisioned  has  the  beam  injected  into  a  very  low  density 


m 
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plasma  In  a  modified  becacron  configuration.  The  question  then  is  how  does 


the  plasma  respond  to  the  beam,  and  more  specifically  how  is  f  related  to 


plasma,  beam  and  system  parameters.  We  treat  the  electron  and  ion  responses 


separately.  Throughout,  we  assume  that  the  background  plasma  has  sufficiently 


small  density  compared  to  the  beam,  that  the  electric  fields  from  the  beam 


dominate  those  from  the  plasma.  Since  the  plasma  is  nearly  at  rest,  there  is 


no  y  cancellation  of  self  electric  fields. 


The  plasma  electrons  react  on  two  time  scales,  the  fast  inertial 


time  scale  and  the  slow  collisional  time  scale.  When  the  beam  enters  the 


plasma,  a  strong  inward  electric  field  is  set  up,  *  -  2irnepip,  where  we  have 


adopted  a  cylindrical  model  for  the  beam  and  p  is  the  radius  (in  the  pololdal 


plane).  Assuming  this  field  is  set  up  slowly  compared  to  an  inverse  cyclotron 


time  (2.5  x  10~1U  sec  for  a  2  KG  field),  the  beam  electrons  respond  by  E  x  B 


drifting  in  the  9  direction  and  drifting  outward  due  to  the  inertial  drift. 


It  is  the  inertial  drift 


„  _  dp  _  -c  dE  .  -c  d  , 

VDi  dt  Boi  dt  Bw  dt  2irabqp 

ce  ce 


which  expels  the  electrons,  (note  q  ■  — | e | )  from  the  beam.  If  an  electron 


starts  out  at  p^  when  n^  ■  t  ■  0,  Eq.  (19)  can  be  integrated  once  to  give 


Equation  (20)  has  an  apparent  divergence,  but  of  course  the  expression  for 


Ep  and  p  are  valid  only  for  p  <  pv,  the  beam  radius.  However,  Eq.  (20)  does 

b  2 

show  that  the  electrons  are  totally  expelled  by  the  beam  if  -■ay  ■  >  1.  For  our 

2w* 

ce 


Thus  the  electron  radius  Increases  exponentially  In  time  with  growth 
2  2 

rate  vw^e/2uce>  Classically 


,  ,  . _  3/2* 

3.4  x  10  T 

e 
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For  temperatures  of  about  1  eV  and  np  ~  10  cm  ,  the  Coulomb  logarithm 
A  -  10,  so  v  ~  3  x  10^.  Thus  the  remaining  electrons  are  expelled  on  a  time 


scale  of  about  20  psec. 

Since  the  electrons  are  forced  away  from  the  beam  on  a  20  psec  time 
scale,  and  even  without  the  beam,  the  electrons  cannot  be  confined  In  a 
toroidal  chamber,  we  expect  the  electrons  to  be  expelled  on  a  time  scale  of 
some  tens  of  microseconds.  This  time  is  long  compared  to  the  time  for  the 
beam  to  center  itself,  but  short  on  the  time  scale  of  the  beam  acceleration. 
Thus,  once  the  beam  begins  to  accelerate,  there  should  be  virtually  no  plasma 
electrons  present. 

For  the  ions  however,  the  story  is  different  because  there  is  a 
strong  attractive  force  between  the  ions  and  the  beam. 

To  continue,  we  consider  the  ion  response.  If  an  ion  is  trapped  near 

.  ■  (2irn,  e^/M)* ^  ~ 


the  center  of  the  beam,  its  oscillation  frequency  is  u 


2  x  10°  sec-1  for  procons  in  our  standard  10  kA  beam.  Since  this  is  ten 
times  larger  than  the  ion  cyclotron  frequency,  the  ions  are  effectively 
unmagnetized.  The  ion  oscillation  time  is  also  much  less  than  the  poloidal 
drift  time  of  the  beam. 

Our  model  is  the  same  as  that  in  a  recent  study  of  the  ion  resonance 
instability,^  namely  that  the  ion  is  initially  at  rest  and  then  when  the  beam 
enters,  the  ion  begins  to  oscillate  due  to  the  electric  field.  Since  ions 
initially  within  the  beam  do  not  leave,  and  other  ions  initially  outside  the 
beam  spend  at  least  part  of  their  oscillation  inside  the  beam,  the  ion  density 
in  the  beam  increases.  We  will  now  estimate  this  increase. 

To  estimate  the  steady  state  ion  density,  we  assume  the  system  is 
cylindrically  symmetric  about  the  beam  center.  Furthermore  all  ions  oscillate 
with  slightly  different  frequencies  so  that  after  several  oscillations,  the 
ions  phase  mix  and  are  distributed  uniformly  along  their  phase  orbits.  Then, 
it  has  been  shown  that  the  ion  distribution  function  in  velocity  space  is 

2n  6(L)  w(H)  M2 

fi(H’L)  "  we  |  34>/3p|  p  ■  p(H)  (23) 

where  the  constants  of  motion  are 


H  -  1/2MV2  +  q  $(p) 


(a) 


(24) 


L  “  mvp 


(b) 


p(H)  is  the  maximum  radius  of  an  ion  with  energy  H,  u(H)  is  the  oscillation 
frequency  of  an  ion  with  energy  H,  $(p)  is  the  electrostatic  potential  of  the 
beam  and  n^  is  the  preformed  plasma  ion  density.  If  the  ion  orbit  is  entirely 


i 

I  within  the  beam 

:  llf1  "  2tvp  (a) 

i 

1  p(H)  -  (2H/Mu2)1/2  (b)  (25) 

'  d)(H)  -  u^.  (c) 

If  the  ion  is  initially  outside  the  beam 

Iffl  “  2imbepb/p  (a) 

(26) 

P(H)  ’  VXP  (b) 

b 

where  H(p^)  ■  q  $(pb).  The  frequency  of  the  ion  outside  the  beam  cannot  be 
computed  in  closed  form.  As  an  approximation  to  it,  use  the  frequency  of  an 
ion  which  rotates  around  the  beam 


w(H)  -  Ul  pb/p(H). 


(27) 


If  p(H)  >  >  r^,  it  is  not  difficult  to  see  that  this  estimate  is  off  by  a 

1  /2 

factor  of  order  [in  p(H)/p b]  • 

The  total  number  of  ions  trapped  in  the  beam  is  then 


Pv 

2ir/  °pdp  {/ 
o  e$ 

I 


H(pJ  H 
b  dH  +  /HdH 


2nlC(H)  io(H)  (M/2 )  1/2 


H(rK)  w||i|p(H)  p(H  -  e*(p))1/2' 


(28) 


II 


The  first  integral,  labeled  I  is  just  irpf  n.  since  it  represents  those  ions 

D  1 

originally  in  the  beam.  The  last  integral  II,  we  only  estimate.  To  do  so, 
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L>> 


a 


p(H)  ai(H)  „  o(H) 

.|2i|  *ui  P, 

Up'p(H)  1  b 


r(H)  ~  d 


(H  -  e*>1/2  ~  fH(d))1/2 


where  d  is  the  maximum  radius  in  the  pololdal  plane.  In  this  case,  the 


integral  II  ~  wn^^  d  (1  +  In  d/p^)  '  so 


N  •  niirpb  fl  +  ~  (1  +  In  d/p^1  ). 


Thus  on  an  ion  oscillation  time  scale,  the  ion  density  in  the  beam  is 
significantly  enhanced.  If  the  beam  is  centered  as  it  is  in  steady  state  so 
d/pfe  ~  5,  the  ion  density  might  be  enhanced  by  nearly  an  order  of  magnitude. 

If  the  beam  is  near  the  liner,  as  it  is  on  injection,  a  reasonable  guess  for  d 
is  the  distance  from  the  liner,  so  d  ~  2pb>  In  this  case  the  ion  density  can 
be  enhanced  by  perhaps  a  factor  of  2. 

Thus  between  the  expulsion  of  plasma  electrons  from  the  beam  region, 
and  partial  trapping  of  plasma  ions  within  the  beam,  the  neutralizing  ion 
density  within  the  beam  should  be  at  least  as  large  as  the  initial  ion  density 
in  the  plasma.  Finally  we  note  that  on  the  time  scale  of  beam  acceleration, 
there  should  be  no  plasma  electrons  in  the  system.  Therefore  the  accelerating 
field  will  not  be  shielded  out  by  any  background  plasma,  and  the  beam  should 
accelerate  as  would  a  vacuum  betatron. 


.  The  Ion  Resonance  and  Ion  Streaming  Instabilit' 


Another  potential  problem  with  the  plasma  assisted  modified  betatron 
is  the  ton  resonance  Instability.  This  Is  a  particular  concern  because  It  Is 
now  established  that  two  other  similar  devices,  HIPAC^  and  SPAC  II ^  were 
disrupted  by  the  ion  resonance  instability.  However  in  both  of  these  devices 
the  beam  nearly  filled  the  chamber,  making  it  particularly  susceptible  to  the 
l  ■  1  instability.  For  the  modified  betatron  with  p.  <  <  a,  the  l  m  1  mode 

D 

should  not  go  unstable  and  the  main  danger  is  an  t  -  2  mode.  This  mode  was 
not  observed  on  HIPAC  or  SPAC  II.  In  the  modified  betatron,  even  if 
parameters  are  right  for  it,  there  is  still  a  good  chance  that  it  will  be 
stabilized  by  the  diffuse  profile.^ 

The  ion  resonance  instability  can  occur  if  the  ion  bounce  frequency 
is  roughly  equal  to  the  1*2  diocotron  mode  frequence.  According  to  Ref.  15, 
this  can  occur  only  if 


For  our  standard  parameters  for  10  kA  beam,  w.  ~  2.4  x  10^,  o>  ~  4  x  10^, 

be  ce 

M/m  ■  1800,  Eq.  (31)  above  reduces  to  y  <  3.  Thus  as  long  as  y  >  3  after  self 
field  diffusion,  the  plasma  background  should  not  give  rise  to  an  ion 
resonance  instability. 

Another  concern  is  the  ion  streaming  instability,  also  discussed  in 
Ref.  15.  In  this  instability  a  cyclotron  mode  on  the  beam  resonates  with  the 
stationary  ions.  The  instability  only  occurs  if  the  parallel  wave  number  is 
in  the  range 


and  the  growth  rate  is  roughly 


n,  a). 

Imu  ~  —  - — - .  (33) 

n  «  u  p,  - 

e  8y2  u 

ce  '  c 

For  our  standard  parameters,  the  growth  time  is  about  20  usee  and  the  range  of 
unstable  k  is  several  percent  and  is  dependent  on  y.  The  idea  then  is  to 
accelerate  the  beam  so  that  y  changes  by  several  percent  in  a  few  growth 
times* 


V.  Conclusions 

We  have  shown  that  a  low  density  background  plasma  can  significantly  aid 
the  modified  betatron  by  allowing  high  current  injection,  avoiding  the 
potential  disruption  on  diamagnetic  to  paramagnetic  transition  and  stabilizing 
the  l  *  2  resistive  wall  instability.  There  are  difficulties  regarding  the 
plasma  production  and  the  ion  resonance  instability.  However  ways  to  overcome 
these  problems  do  exist  and  appear  to  be  feasible. 
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DIAMAGNETIC  BEAM 


a)  The  dependence  of  on  radius  for  a  diamagnetic  beam.  At 


the  reference  orbit,  is  a  relative  minimum. 

b)  Same  for  a  paramagnetic  beam.  At  the  reference  orbit,  Pa  is 

9 


a  relative  maximum. 


2: 


In  this  appendix  we  discuss  some  of  the  recent  simulations  of  a  beam  in 
the  modified  betatron  as  it  is  accelerated  through  the  diamagnetic  to 
paramagnetic  transition.  In  Ref.  18,  a  10  kA  beam  was  accelerated  from  3  MeV 
to  10  MeV  by  linearly  ramping  up  the  vertical  field  over  a  time  of  about  500 
nsec.  If  the  beam  has  no  spread  in  longitudinal  energy,  it  seems  to  be  able 
to  survive  this  transition  with  about  a  30Z  Increase  In  emlttance  and  beam 
radius.  On  the  other  hand,  with  a  1Z  spread  in  longitudinal  energy,  the 
emlttance  increased  by  about  an  order  of  magnitude  and  the  beam  radius  also 
increased  significantly.  In  addition,  for  the  latter  case,  the  phase  space 
plots  clearly  show  particles  scraped  off  the  beam  edge  and  drifting  out  along 
a  separatrix. 

The  simulation,  at  least  the  first  (constant  longitudinal  energy)  case, 
is  considerably  more  optimistic  than  what  the  fluid  theory  of  Refs.  4-6  would 
predict.  The  difference  in  the  two  results  probably  has  its  origin  in  the 
difference  in  time  scales.  In  both  the  fluid  theory  and  the  simulations,  P 
is  an  exact  constant  of  the  motion.  However  the  fluid  theory  has  one 
additional  constant  of  motion,  the  toroidal  flux  through  the  drift  orbit.  The 
significant  point  is  that  this  flux  is  not  an  exact  constant,  but  an  adiabatic 
invariant.  Therefore  for  it  to  be  conserved,  the  changes  in  the  orbit  have  to 
be  very  slow  compared  to  the  poloidal  drift  orbit  time.  If  the  time  for 
changes  are  comparable  to  a  drift  orbit  time,  the  whole  concept  of  an 
adiabatic  fluid  evolution  of  these  drift  orbits  breaks  down. 

For  the  modified  betatron,  the  drift  orbit  time  is  about  200  nsec  or 
more.  In  the  simulations,  the  energy  went  from  3  MeV  to  10  MeV  in  500  nsec, 
so  it  is  extremely  unlikely  that  an  adiabatic  theory  would  be  valid.  Since 
there  is  only  one  constant  of  motion,  PQ  ,  the  transverse  orbits  would  most 
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likely  fill  up  the  beam  region  ergodically.  However  in  the  modified  betatron 
experiment,  the  acceleration  time  is  milliseconds.  Thus,  at  least  up  to 
transition,  and  shortly  after  transition,  the  adiabatic  theory  should  be 
valid.  Of  course  adiabatic  theory  cannot  be  valid  at  transition  itself,  as 
pointed  out  in  Refs.  4-6,  and  it  may  be  that  even  for  the  slow  time  scale 
acceleration,  the  beam  may  pull  itself  through  transition  much  more  quickly 
than  the  time  scale  indicated  by  change  in  vertical  field.  Results  from  the 
fast  time  scale  simulations  indicate  that  this  is  the  case.  However  at  least 
during  the  early  phase  of  transition,  the  adiabatic  theory  does  indicate  that 
the  outer  part  of  the  beam  is  scraped  off  onto  a  separatrix  whose  size, 
depending  on  the  parameters,  can  be  comparable  to  the  liner  radius. 
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